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Abstract. Using results on the mass gap in the sine-Gordon model combined with the exact amplitudes 
in the bosonized representation of the Heisenberg spin-1/2 chain and one-loop renormalization group, we 
derive a quantitative expression for the gap in a dimerized spin-1/2 chain. This expression is shown to be 
in good agreement with recent numerical estimates when a marginally irrelevant perturbation is taken into 
account. 

PACS. 75.10.Pq spin chain models 



Low dimensional antiferromagnets have been the sub- 
ject of intense scrutiny both theoretical and experimental 
for the last twenty years. The simplest model, the spin- 

' 1/2 Heisenberg antiferromagnet, is integrablepQ| and can 
be mapped onto a continuum field theory which al- 
lows the full determination of its zero temperature critical 
behavior. The presence of a marginally irrelevant operator 
in the continuum theory induces logarithmic corrections 
to the critical scaling^. The corrections to scaling of the 
correlation functionsjSHEllZj, NMR relaxation rates [HUH], 
and susceptibilities [lOllll| in this model have been inves- 
tigated in details. Further, when the Heisenberg spin-1/2 
chain model is perturbed by a relevant operator such as 
an alternation of the exchange coupling, the marginal op- 
erator gives rise to a logarithmic correction to the power 
law dependence of the gap on the perturbation 5 . Such 
logarithmic corrections to scaling in the gap in the con- 
text of two dimensional statistical mechanics of the four 
state Potts model whose transfer matrix is related to the 
Hamiltonian of the alternating Heisenberg chain [T2Tll3lll4l 

■ UK] , In m particular, it was shown that the dependence 
of the gap A on the dimerization S, was changed from the 
form A ~ <5 2 / 3 [WP7] to a the form A ~ 8 2 ^ /\ ln^l 1 / 2 . 
Such logarithmic behavior was confirmed by numerical 
calculations in |18II19II20| . Alternatively, the dependence 
of the gap on the dimerization can be described by an 
effective power law form with an exponent that deviates 
from 2/3 |21U22| . For a dimerization not too small, it is 
found that the resulting effective exponent is close to 2/3 
|21|. Further, by considering a Heisenberg chain with an 
additional next-nearest neighbor coupling finely tuned to 
cancel the marginal operator, a pure power law with ex- 
ponent 2/3 can be obtained obtained for the gap |23| . 
Recently, the logarithmic corrections were investigated in 
greater details using the DMRG[2~4"). The data for the gap 



could be fitted to the form: 



A = al'l 



£2/3 



9ap (In So/5) 1 / 2 '' 



(1) 



with a gap = 19.4 and <5o = 115 or alternatively by the 
power law form A = 1.94<5 073 . A difficulty that arises 
when comparing the predictions of the Renormalization 
Group approach [51ll2lll3lll4lll5j with the numerical results 
is that the former approach can only predict the expo- 
nents, and not the non-universal prefactors. However, ex- 
act results for the sine-Gordon model combined with re- 
cent progress [25I27S] on the bosonization treatment of the 
Heisenberg spin 1/2 chain using integrability make it pos- 
sible to overcome these two difficulties and obtain the pref- 
actor in the expression of the gap at least in the absence 
of logarithmic corrections. Assuming that the gap varies 
continuously as the marginally irrelevant is turned on, it is 
then possible to obtain an expression of the gap as a func- 
tion of the marginally irrelevant operator, with no further 
unknowns. Fitting the data of |2U then allows the determi- 
nation of the order of magnitude of the marginally irrele- 
vant interaction. The obtained value can then be checked 
against the one obtained in A similar approach has 
been used previously in _ 2"7] to estimate the gaps induced 
by a staggered field in an anisotropic spin 1/2 chain. 
The Hamiltonian of the dimerized spin-1/2 chain reads: 



H= J^(l + (-)"5)S n -S n+1 



(2) 



For 5 = 0, this Hamiltonian reduces to the one of the uni- 
form antiferromagnetic Heisenberg chain (J > 0) the low 
energy properties of which are described by the following 
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H = / 



252 



continuum Hamiltonian j4~lll4| : 

uK{nII) 2 

where [<^(x), JI(x')] = z<5(x — a;'), it 

The latter conditions ensures SU{2) symmetry, and 



dx 
2^ 



(27ra) 2 



Heaviside step function. Finally,n 2 (a; + a)(Ji + Jl){x) 
and n z (x)( Jl + Jt )(x + a) are respectively obtained from 
dxcos VSxfcffi^ 8 ' ® an< ^ ® ^ e tw0 f° uowm g short distance expan- 



= \ Ja and X = 



1 - 



for 52 < the operators (nil) 2 — (d x <f>) 2 and cos are 
both marginally irrelevant resulting a gapless fixed point. 
The exact bare value of gi has been estimated in 

The spin operators can be expressed as a function of 
as: 



17, 



— = (J+ + J_)(no) + (-)' T/a n( M ), 



J+(as) = (J? + iJ r y )(x) = — e 

27ra 



-i^2(fi-r<f)(x) 



i+(x) = 



rll - d x (, 



2ttV2 
A 



(x) = — $m^/24> a {x), 



(4) 

mvi, (5) 

(6) 

(7) 
(8) 



A = 



(?) 



1/4 



(9) 



In order to determine the bosonized Hamiltonian of the 
dimerized spin-1/2 chain (J2J), all we need is a bosonized 
expression of the dimerization operator ^2 n (— )™S„ • S„+i. 
Using J3J , the corresponding expression is easily extracted 
from: 



is 

a 



S n+ i = (uniform) + 

{-) n a [-(J+ + J_)(no) • n((n + l)a) 
+n(no)-(J + + J_)((n + l)o)] (10) 



The bosonized expression of the dimerization operator, is 
thus obtained from the the short distance expansion of 
Jfl,z, and n. Using Eqs. JSJ) and with Glauber identi- 
ties, one finds the following expressions: 



n ± (x)(J+ + J-) T (x 
(,/+ + J-) ± {x)7i f {x 



+ a) = 



A 



(tto) 2 

a) = — — T7 
{Tray 



os \/2(j){x) + 
■ cos x - 



.(.11a) 
(lib) 



The change of sign is a consequence of the application 
of the Glauber identity taking into account the commu- 
tation relation [cj>(x),9(x')] = mY(x' — x), Y being the 



1 



2tt\/2 
1 

2W2 



d x (j)(x + a) sin V2(f>(x) — — !— cos V2(f)(x) +(.12a) 

2ira 



d x <p(x) sin y2(f)(x 



a) = --J—cosV^^x) ft2b) 
27ra 



which can be derived by normal ordering the product of 
the two operators |29U3U| . A sketch of the derivation is 
given in the appendix. It is easily seen that Eqs. (|llf> - (|12f) 
are compatible with spin rotational invariance. Combining 
the expressions l|ll|l and 1)12(1 in l)10[l. and using the value 
of A in Eq. wc finally obtain that: 



1 



S„ • S„+i = uniform + (-) r 



(f) 



1/4 



cos 



^(13) 



Therefore, the continuum Hamiltonian describing the dimer- 
ized spin 1/2 chain at low energy reads: 



H = / =L 



where a is a lattice spacing, jj-j- represent Majorana fermion 
operators that can be omitted in some cases (see e. g. Ref. 
|2*8) for a discussion of this point), and 9 is defined by 
6{x) = 7T f*^ dx'n(x'). The constant A is a non-universal 
parameter that depends on the lattice model being con- 
sidered. Recently, this parameter has been determined in 
the case of the isotropic Heisenberg spin-1/2 chain 2 5Tl26| . 
and it was found that: 



dx 

2ff 2 
' (27ra) 2 



dx cos 



2gi 



(27T0) ; 



■ cos 



V2<j) 



(14) 



Note that in 1)14(1 , the sign of g± does not matter as it can 
always been rendered positive by the shift cj) — > cfr + ir/v2. 
In 1(14(1 . we have: 



91 



= 6J 



(f 



1/4 



5 a. 



(15) 



52 is a marginally irrelevant field 
0. Let us assume for a moment that 



As we noted before 
which flows to if g\ 
we can neglect completely the presence of this marginally 
irrelevant operator and take K = 1, 32 = in 1(14(1 . Then, 
the Hamiltonian 1(14(1 becomes the sine-Gordon model. 
This model is integrable, and the expression of the gap 
can be found in or in Eq. (12). In the notations 
of H21, P 2 = K/A = 1/4, and /x = 3/ 7 r 3 (7r/2) 1 / 4 ( 5 (where 
we have used the fact that the velocity u — \Ja). The 
dimcnsionless gap M is then given by: 



M 



2 r(i/6) 

7^(2/3) 



r(3/4) 3 ^7ry/4 

r(i/4)^ 



(5) 



and the energy gap is given by A 

1.723<5 2/3 



2/3 



M i.e. 



A 
1 



1.097S 2 (%6) 



(17) 



We note that the formula (|TBjl has already been applied to 
calculate the gap of the dimerized spin 1/2 chain in |33j . 
but the value of A, Eq. 10 was not known. The formula 
(117(1 is in reasonable agreement with the result quoted in 
[2*T] who reported that A/ J = 1.5<5 ' 65 as two expression 
differ at most by 6% for 0.01 < 6 < 0.1. Comparing our 
expression flJJ to the one of Ref. EH, A/ J = 1.94<5 a73 , 
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Fig. 1. Comparison of the gap obtained in Ref . |24) (solid 
line) with the expressions (|17|l without logarithmic correc- 
tions(dashed line) and the expression 1291 including logarith- 
mic corrections with 2/2(0) = —0.22 i.e. a = 4.499 and So = 148 
(dash dotted line). Equation 117H gives a result nearly indis- 
tinguishable from the result of Ref. |24). 
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Fig. 2. Comparison of the expression of ground state energy 
of Ref. [24 ] with the expression derived in this paper without 
logarithmic corrections 118H (dashed line) and with logarithmic 
corrections (gOl for y 2 (0) = -0.22 i.e. (3 = 1.86 and S = 148 
(dash-dotted line) . Including logarithmic corrections leads to a 
better agreement for small dimerization. However, deviations 
are still significant in contrast with the case of the gap. 



we find that they are in agreement within a 10% relative 
error when 5 > 0.03 as represented on Fig. For lower 
values of 5, the two results deviate sensibly. As we shall 
see, this is the result of the logarithmic corrections. 

In , the expression of the ground state energy was 
also given in dimensionless units in Eq. (14). Using this 
expression, we obtain for the ground state energy: 



E __7T 

J ~ 2 J 



M 2 IT 

tan — 

4 6 



-0.2728<5 4/3 



(18) 



This expression is compared to the one quoted in |24| , 
E / J = -0.398 1 - 45 on Fig. El The two energy formulas 
are in better agreement than the gap formulas at low 
dimerization. Till now, we have totally neglected the pres- 
ence of the marginally irrelevant operator cos %/80- As we 
shall now see, the corrections to scaling [THIS] induced by 
this operator in the gap formula, are responsible for the 
discrepancies between the numerical and the analytical 
results. The renormalization group equations associated 
with the Hamiltonian (|14|l read|13|: 



d 
Jl 



1 

K 
dyi 
dl 

dm 
dl 



1 2 1 2 

gZ/i + 2^2, 



K 



= I 2- — + 2/2 ] yi, 



(2 - 2K)y 2 



(19) 
(20) 
(21) 



where we have introduced yi = gi/(iru). For y\ = 0, the 
SU(2) symmetric flow is recovered for K = 1 — 2/2/2. 
Then, the equations l|19fl reduce to the single Kosterlitz- 
Thouless|3H| dy 2 /dl = y%. We see that for y 2 < 0, this 
equation flows to the fixed point y 2 — 0, with the follow- 



ing flow equation: 



2/2(0 



V2 (0) 



i-y 2 (o)z 



(22) 



Let us now assume |T3"] that we have turned on a very small 
2/1. Using the initial conditions with SU(2) symmetry, we 
can easily show that the RG equations reduce to: 



dy 2 
dl 

dyi I ■-> 

dl 



2 1 1 2 
2/2 + » 



2/2 yi. 



(23) 
(24) 



If we assume that yi(0) <C 2/2(0), we can assume that in 
(|23fl . we can take y% = 0, so that the flow of y 2 is given 
by H22|l. Then, the equation Ij24(l is trivially integrated, 
leading to: 



2/i(0-2/i(0) 



(1 + |2/ 2 (0)|0 3/4 



(25) 



This equation should break down for a scale lo such that 
2/i ('0) ~ 2/2 (lo)- One has: 



e io (l + 2/ 2 (0)?o) 



1/6 _ 



2/3 



2/i(Q) 2 /3 



(26) 



For I > Iq the contribution of y 2 to the renormalization of 
yi being negligible, yi(l) = e 3 / 2 ^~ io) 2/i(/ ). The scale /* 
at which 2/1 (0 ~ 1 ^ s thus given by: 



(l + |2/ 2 (0)|Zo) 2 / 3 
|2/ 2 (0)| 2 / 3 



l2/i(0)| 



2/3 



(1 + |2/ 2 (0)|Z ) 1 / 2 



(27) 
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An approximate form of lo can be obtained by iterating 
(J2EI leading to: 



l!/i(0)| 



2/3 



1 



2/2(0)| In 



1/2(0) 

yi(o) 



1/2 



(28) 



Since the scaling of the gap is A ~ e and the scaling 
of the ground state energy is E ~ e _2/ , for 2/1 — >■ 0, these 
formulas are in agreement with the scaling predicted in 0. 

We now make an important assumption. We assume 
that in the formulas l|17|) (|18|l . we can replace M with 
Ce~ 1 ', I* being given by (|28|l . with C being chosen in such 
a way that for y% = 0, the resulting M agrees with (116(1 . 
This is a uncontrolled approximation as the model (|14fl is 
non-integrable but it is partially justified by the fact that 
the energy and the correlation length evolve continuously 
as a function of the parameter y% in the vicinity of the 
integrable point. With this assumption, and noting that 
the definition of j/i and Eq. I|15fl imply y\ = 1.3612(5, we 
obtain that C = 0.8932. 

We are lead to the following expressions of the gap: 



A 
7 



1.7235 2 / 3 



1 



y 2 (0)| In 



2/2(0) 



1.36125 



1/2 : 



(29) 



and of the ground state energy difference: 
E 0.27286 4/3 



J 



1 



2/2(0)| In 



3/2(0) 



1.36125 



(30) 



The author thanks R. Chitra, R. Citro, T. Giamarchi and P. 
Lecheminant for discussions. 



A derivation of the short distance expansion 

The short distance expansion (|12fl has been derived in [291 
I30j using a fermionic representation. In this appendix, we 
give an alternative derivation using the bosonic represen- 
tation. We start from the equation: 

: d x <j>{x') :: V{<j>{x)) :=: d x <P(x')V(<P(x)) : + : ^(</>(x)) : d x {<j>(x')ct>{x) 



which is easily obtained by expanding V(cj>) as a power 
series and applying Wick's theorem|3"5|. In the case of the 
massless free boson, we have: 



These equations are trivially reduced to the Eqs. (7) and 
(8) in [21], by expressing a, a gap , So as a function of 2/2(0). 
A good fit to the power-law expression of the gap quoted 
in [21] is obtained using Eq. J2U with 2/2(0) = -0.22 (see 
Fig.f^). The corresponding fit for the ground state energy 
using Eq. (|30|l with the same value of 2/2(0) = —0.22 is bet- 
ter than the fit obtained without logarithmic corrections 
(see Fig. |2J , however it is not as good as the fit obtained 
for the gap, especially for 6 > 0.01. This could be due 
to the non-singular part of the free-energy which is not 
taken into account in the Renormalization Group calcu- 
lation. The value of 2/2(0) = —0.22 compares reasonably 
well with the one quoted in |51l35| 2/0 = —0.25. 

To summarize, we have shown that the results of Ref. 24 
could be recovered from a bosonization approach includ- 
ing the appropriate operator renormalizations, and using 
exact results for the sine-Gordon model combined with a 
one-loop RG. The amplitude of the marginally relevant 
operator was found to be in reasonable agreement with 
an independent estimate coming from logarithmic correc- 
tions to the dependence of the gaps in a spin-1/2 chain. 
Given the relatively large value of the coupling constant, 
we are at the limit of applicability of the one-loop RG. 
Better agreement might be obtained by going beyond the 
one-loop approximation 3(5 • The present approach does 
not depend crucially on integrability as it is also possible 
to determine the parameter A in |TJJ for a non-integrable 
model via numerical calculations |27U37| . 



1 



x — x 



(d>(x')cf>(x) - d> 2 ) = -In 
which leads to the expansion: 
: d x( j){x') :: V{(j){x)) :=: d x dp{x')V(dp{x)) : + 



2(x' - x) 



(32) 



5H0(z)O33) 



Applying this formula in the case of V(<p) — sin \[2<$> leads 
to Eqs. O. 
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